Continuous quantum feedback of coherent oscillations in a solid-state qubit 
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We have analyzed theoretically the operation of the Bayesian quantum feedback of a solid-state 
qubit, designed to maintain perfect coherent oscillations in the qubit for arbitrarily long time. 
In particular, we have studied the feedback efficiency in presence of dephasing environment and 
detector nonideality. Also, we have analyzed the effect of qubit parameter deviations and studied 
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I. INTRODUCTION 

Continuous quantum feedback in optics and atomic 
physics has been studied theoret icall'v^ ^'^i for more 
than a decade 1 (see also Refs. kl7lglfllloi ) and has been 
recently demonstrated experimentally^ In contrast, con- 
tinuous quantum feedback in solid-state mesoscopics is a 
relatively new subject The use of quantum 
feedback to maintain coherent (Rabi) oscillations in a 
qubit for arbitrarily lon g t ime has been proposed and 
analyzed in Refs. Il3 an dll3t a simplified experiment has 
been proposed in Ref. Il4l Cooling of a nanoresonator 
by quantum feedback has been proposed and analyzed in 
Ref. The use of quantum feedback for the nanores- 
onator squeezing has been studied in Ref. [la- 
Feedback control of a quantum system requires con- 
tinuous monitoring of its evolution (in ideal case the 
wavefunction should be monitored), which is the main 
non-trivial feature of quantum feedback. Obviously, the 
operation of quantum feedback cannot be analyzed us- 
ing the "orthodox" approach of instantaneous collapse, 17 
which is not suitable for continuous quantum measure- 
ment. Also, the ensemble-averaged ("conventional") 
approach^ to continuous quantum measurement is not 
suitable since it cannot describe random evolution of a 
single quantum system. Therefore, analysis of quantum 
feedback requires a special theory capable of describing 
continuous measurement of a single quantum system. 

The development of such theories has started long 
agoA2i22i2i and has attracted most of attention in quan- 
tum opticsi*22i2^ (in spite of similar underlying princi- 
ples, the theories may differ significantly in formalism 
and area of application) . For solid-state qubits such the- 
ory ( "Bayesian" formalism) has been developed relatively 
recentlji2£i2ii (for review see Ref. l26j) . The equivalence of 
the Bayesian formalism to the quantum trajectory ap- 
proach translated^LSSiS from quantum optics has been 
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shown in Ref. |27[ 

In simple words, the Bayesian formalism takes into 
account the information contained in the noisy output 
of the solid-state detector measuring the qubit, so that 
the corresponding quantum back-action onto qubit evo- 
lution is described explicitly. In classical probability the- 
ory the way to deal with an incomplete information is 
via the Bayes formula^ it can be shown24i2&i2& that a 
somewhat similar procedure should be used for evolution 
of the qubit density matrix due to continuous measure- 
ment, that explains why the formalism is called Bayesian. 
The Bayesian formalism shows that an ideal (quantum- 
limited) detector can monitor precisely the random evo- 
lution of the qubit wavefunction in the course of measure- 
ment; and if the measurement starts with a mixed state, 
the qubit density matrix is gradually purified, eventu- 
ally approaching a pure state. The quantum point con- 
tact (QPC) is an example of (theoretically) ideal detec- 
tor. When the detector does not have 100% quantum 
efficiency [as in the case of a single-electron transistor 
(SET)] , there is an extra dephasing term in the evolution 
equation, so that the qubit purification due to gradually 
acquired information competes with the decoherence due 
to detector nonideality. 

The possibility to monitor the random quantum evo- 
lution of the qubit in the process of measurement natu- 
rally allows us to arrange a feedback loop which keeps the 
qubit evolution close to a desired "trajectory" . Of course, 
the measurement process disturbs the qubit evolution; 
however, the detector output contains enough informa- 
tion to monitor and undo the effect of this disturbance. 
It is important that the deviations from the desired tra- 
jectory due to interaction with decohering environment 
are efficiently suppressed by the feedback loop, that can 
be useful, for example, in a quantum computer. The feed- 
back loop considered in Refs . l25l andll^ has been designed 
to maintain the coherent oscillations in the qubit for arbi- 
trarily long time by comparing the oscillation phase with 
the desired value and keeping the phase difference close 
to zero (the amplitude of oscillations is equal to unity in 
the case of ideal detector and energy-symmetric qubit). 
It has been shown that the fidelity of such feedback loop 
can be arbitrarily close to 100%, while it decreases in 
the case of a nonideal detector and/or significant inter- 
action with environment as well as in the case of finite 
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FIG. 1: Schematic of the quantum feedback loop maintaining 
the quantum oscillations in a qubit. The qubit oscillations 
affect the current I(t) through a weakly coupled detector; 
this signal is translated by the "processor" into continuously 
monitored value p m (t) of the qubit density matrix. Next, by 
comparing p m (t) with the desired oscillating state p (t) 7 a 
certain algorithm ("controller") produces the feedback signal 
applied to an "actuator" which changes the qubit tunneling 
amplitude H+AHfb, in order to reduce the difference between 
p m and p d . 



bandwidth of the line carrying the signal from detector. 

The present paper is a more detailed analysis of the 
operation of the feedback loop proposed in Refs. and 
H2T In particular, we study the feedback loop operation 
in presence of extra dephasing due to environment and 
nonideal detector, analyze the effect of qubit parameter 
deviation, and consider the feedback of a qubit with en- 
ergy asymmetry. In the next Section we describe the 
model, in Section III we consider the feedback operation 
in the ideal case, Section IV is devoted to the effects of 
nonideal detector and extra dephasing, in Section V we 
analyze the worsening of feedback efficiency in the case 
of qubit parameter deviations, in Section VI we study 
the feedback of an energy- asymmetric qubit, and Section 
VII is a conclusion. 



II. MODEL 

Let us consider the quantum feedback loop shown in 
Fig. ^ which controls the qubit characterized by the 
Hamiltonian 



H 



qb 



£ fb 1 t 



cjci) + H fb (c[c 2 + 4 c l)) (!) 



where c] 2 an d c i,2 are creation and annihilation opera- 
tors corresponding to two "localized" states of the qubit, 
representing the "measurement basis" . The qubit energy 
asymmetry En and tunneling amplitude can both 
be controlled by the feedback loop: 



H fb =H + AH 



Sfb 



Ae 



(2) 



however, in this paper we assume Aef b = 0, so that only 
tunneling is controlled. Intrinsic frequency of coherent 
oscillations in the qubit (without interaction and feed- 
back) is f2 = V '477 2 + e 2 /h; we also call it Rabi frequency, 
not implying presence of microwave radiation (despite 



this terminology differs from the initial meaning of Rabi 
oscillations, it is conventionally used nowadays). 

For definiteness we consider a "charge" qubit continu- 
ously measured by QPC or SET, so that the measure- 
ment setup is similar to what has been studied thco- 

Taking into account the quantum back-action due to 
measurement, the evolution of the qubit density ma- 
trix p is described by the Bayesian equationa2^2^i2& (in 
Stratonovich form) 

77 2AZ 
pu = -2 Imp^ + P11P22 -PT~ - Mi ( 3 ) 
n bi 

■ £ fb ■ Hfb 1 \ 

Pi2 = i-r- P12 + 1 — =- (pn ^ P22) 
n n 

-{pu - ^22)^- [I(t) - Jo] pn - 7P12 , (4) 

where I(t) is the noisy detector current (output signal), 
Si is the spectral density of current noise, AI = 7i — I 2 is 
the difference between two average currents 7j and 7 2 cor- 
responding to the two qubit states, and 7 = (ii + 
The dephasing rate 7 = 7^ + j env has the contribution 
7d due to detector nonideality, — — l)(A/) 2 /4Sj 
(here 77 < 1 is the quantum efficienc y 24 ' 2 ^' 2 ^'^^ 4 '^ ) and 
contribution j en v due to interaction with extra environ- 
ment. As always, pu + P22 = 1 and p%i — p*2- Equations 
(J3J (@J imply weak detector response |A7| <C 7 , quasi- 
continuous current, and large detector voltage compared 
to the qubit energy. The current 



i(t) = i a + (Ai/2)( P11 -p 22 )+m 



(•5) 



has the pure noise contribution £(t) with frequency- 
independent spectral density Sj. Notice that averag- 
ing of Eqs. ©~0 over leads to the standard 
ensemble-averaged equations^ with ensemble dephasing 
rate T = (A7) 2 /ASi + 7. We characterize coupling be- 
tween qubit and detector by the dimensionless constant 
C = h(AI) 2 /SiH (we assume^ 77 > 0) and concentrate 
on the case of weak coupling C < 1 (notice that C = 1 
can still be considered a weak coupling since the qual- 
ity factor of oscillations in presence of measurements^! is 
877/C for e = 0). 

In this paper we consider the "Bayesian" feedback^ 
which requires a "processor" solving Eqs. in real 

time - see Fig. ^ (other possibilities are, for example, 
"direct" feedback briefly mentioned in Ref. and "sim- 
ple" feedback via quadrature components analyzed in 
Ref. IT^l . In this paper we neglect the effect of finite 
bandwidthi^iii^ of the line carrying the detector signal, 
and we also neglect the signal delay in the feedback loop. 
As a result, in most of the paper we assume that the 
monitored value p m (t) of the qubit density matrix coin- 
cides with the actual value pit). Only in Section V we 
consider p m different from p because of the deviation of 
the qubit parameters 77 and e from the values assumed in 
the "processor" (finite signal bandwidth would also lead 
to difference between p and p m ). 
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For the feedback control the monitored qubit evolution 
is compared with the desired evolution (Fig. , and the 
difference signal is used to control the qubit parameters 
in order to decrease the difference. Actually, various al- 
gorithms ("controllers") are possible for this purpose; in 
this paper we will consider linear control (see below). We 
study the feedback loop, which goal is to maintain per- 
fect coherent oscillations in the qubit for arbitrarily long 
time, and (except for Section VI) the desired evolution is 



d 

Pll 



i + cos riot 



d sin Clot 

Pl2 = * 7, . 



(0) 



with frequency = 2H/U corresponding to e = 0. Ex- 
cept for Section VI, we assume the following feedback 
control: 

AH fb = -FHA4> m , (7) 

Acj> m = 0m (t) - fl Q t (mod27r), (8) 
<Mi) = arctan[2Imp™/(Pn - P&)] 

+(7r/2)[l- Sg n(pE-pg)], (9) 

where 4> m is the monitored value of the phase, phase dif- 
ference A<fi m is defined as |A0 m | < tt, and F is a di- 
mensionless feedback factor [the second term in Eq. © 
provides proper phase continuity on 2tt circle] . The con- 
troller Q is supposed to decrease the phase difference 
(negative feedback): if phase <p m (t) is ahead of the de- 
sired value, then AH ft, is negative, that slows down the 
qubit oscillations and decreases the phase shift; if 4> m (t) 
is behind the desired value, the oscillation frequency in- 
creases to catch up. 

We will characterize the feedback efficiency by the 
"synchronization degree" D defined as averaged over 
time scalar product of two Bloch vectors corresponding 
to the desired and actual states of the qubit. An equiva- 
lent definition is 



D = 2(Trpp d ) - 1, 



(10) 



where (..) denotes averaging over time. Perfect feedback 
operation corresponds to D = 1 (notice that p d is a pure 
state). Feedback efficiency D can be easily translated 
into average fidelity as (D + l)/2 or \J {D + l)/2, de- 
pending on the definition of fidelitjiifi^ (translation for- 
mula would be slightly longer if neither p nor p d are pure 
states). We prefer to use D instead of fidelity because 
D = in absence of feedback when p and p d are com- 
pletely uncorrelated, while fidelity is non-zero. 



III. IDEAL CASE 

Let us start analysis with the basic ideal case of r) = 1 
(quantum- limited detector, e.g. QPC), absence of extra 
environment (j env — 0), and symmetric qubit (e = 0). 
The analytical results for this case have been presented 
in Ref. 12| here we discuss the derivation in more detail. 



Since rj — 1 and "f en v = 0, so that there is no dephas- 
ing term in Eq. Q, the qubit density matrix p becomes 
pure in the process of measurement Because of the 
energy symmetry, efb = £ = 0, the real part of pi2 even- 
tually becomes zero. This happens because the product 
(pii — Pti){J — Iq) affecting the evolution of Repi2 in 
Eq. (@J is on average positive. Therefore, after a tran- 
sient period the evolution of the density matrix p can be 
parameterized as 



pn = (1 + cos0)/2, pi2 = z(sin<£)/2 



(11) 



with only one parameter cj>(t). We have also checked this 
fact numerically. Notice that since the qubit is monitored 
exactly, p = p m , the phase <f> coincides (modulo 2tt) with 
the monitored phase </> m defined by Eq. JJjJ. 

The evolution equation for phase <f> can be easily de- 
rived from Eq. as 



2H fb /h- (A7/5j)(/-/ )sin< 



(12) 



so the phase difference A(f> 
with A<f> m ) evolves as 



flat (which coincides 



j t A<f> = -sin^^ COS0 + A -FQ A<f>. (13) 

(All equations are in the Stratonovich form, so we use 
usual rules for derivatives 41) Notice that because of our 
definition \A<p\ < tt, the phase difference jumps by ±2tt 
at the borders of ±7r interval. 

For weak coupling (C/8 <C 1) the qubit oscillations 
are only slightly perturbed by measurement and corre- 
sponding phase diffusion is relatively slow. Assuming 
that the feedback control is also slow on the timescale 
of oscillations (\AHfb\ -C H), we can average Eq. 113|) 
over relatively fast oscillations. Then the first term in 
parentheses is averaged to zero and averaging of the term 
— (sm<j>)(AI / Si)£(t) leads to the effective noise with 
spectral density = (A/) 2 /2Si, so that the remaining 
slow evolution of phase difference is 



dt 



A(j) = £- Fn A(j). 



(14) 



To find the feedback efficiency D = (cosA0) analyt- 
ically, let us also assume that feedback performance is 
good enough to keep the phase difference A</> well inside 
the ±7r interval, so that the phase slips (jumps of A<p by 
±27r) occur sufficiently rare. In this case we can consider 
Eq. (|14fl on the infinite interval of Acj). The corresponding 
Fokker-Planck-Kolmogorov equation for the probability 
density <j(A<p) 



da 



d 



dAcj) 



(aFn Acj>) 



1 d*{S~ e ) 
4 <9(A</>) 2 



(15) 



has the Gaussian stationary solution a s t(A(f>) 
(27rV)" 1/2 exp[-(A(/ ) ) 2 /2V] with variance V 
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FIG. 2: Solid lines: quantum feedback efficiency D as a func- 
tion of the feedback strength F for five different values of 
the coupling C and ideal operation conditions (see text). The 
curves for C < 3 practically coincide with each other. Dashed 
line shows the analytical result l|16^ . Inset shows the same 
curves for larger range of F/C. 



S|/4Ffi = C/16F. Therefore, (cosA0) = exp(-V/2), 
and so the feedback efficiency ia~ 

D = cxp(-C/32F) (16) 

in the case of weak coupling and sufficiently efficient feed- 
back (C < 1, D > 1/2). 

Figure shows comparison between the analytical re- 
sult (|16J) and numerical results for D as a function of feed- 
back factor F (scaled by coupling C). Numerical results 
have been obtained by direct simulation of the Bayesian 
equations ©-(0J using the Monte Carlo method^*2^ for 
five values of coupling: C = 10, 3, 1, 0.3, and 0.1. One can 
see that for weak coupling C < 1 the analytics works very 
well when the feedback is sufficiently efficient, D > 0.5. 
Another important observation is that with the feed- 
back factor F normalized by coupling C, the curves for 
C < 1 are practically indistinguishable from each other, 
the curve for C — 3 goes a little higher but still within the 
line thickness, and only the curve for C — 10 is noticeably 
different. Therefore, as expected, the weak-coupling limit 
is practically reached starting with C < 1. This makes 
unnecessary to analyze numerically the case of very small 
coupling C -C 1, which requires much longer simulation 
time than the case of moderately small coupling. 

Notice that \AHfb\/H < ttF, and F scales with cou- 
pling C. Therefore, in the experimentally realistic case 
C « 1 a typical amount of the parameter change due 
to feedback is small, \AHfb\ <C H. [Hence, we should 
not worry about unnatural assumption of using control 
equation J7J even when Hfb becomes negative.] 

The feedback efficiency D is directly related^ to the 
average in-phase quadrature component of the detector 
current, cosfi *) = (A7/4)[Z»+(cos(2O o *+A0))], so 
that in the case of practically harmonic oscillations D — 
(4/ AI)(I(t) cosilot). Positive in-phase quadrature is one 
of easy ways to verify the quantum feedback operation 



experimentally. 

Besides analyzing feedback efficiency D, let us also cal- 
culate the qubit correlation function K z (t) — (z(t + 
t) z(t)) where z — pn — p22- In the case of practi- 
cally harmonic (weakly disturbed) oscillations, the corre- 
lation function K z (t) = (cos[0(i + r)] cos[0(i)]) is equal 
to (cos[£V + <J<£(r)]) /2 where 8<P{t) = A<j)(t + r) - A<f>{t) 
is the phase deviation during time r. Since in our case 
(sin(5(/)(r)) = because of the symmetry of Eq. I|14fl . the 
correlation function is reduced to 



cosfV 
K A T ) = — o — (cos^(r)). 



(17) 



We can find (cos50(r)) using exact solution of the 
Fokker-Planck-Kolmogorov equation l|15|l with initial 
condition <r(A0, 0) = S(A<j) - A</> ): 



ct(A0,t|A0 o ) 



exp[-(A</> - A^e- MlaT f/2V{T)] 



V(t) - (S S ~/4FQ ) (1 



-2Ffl T 



,(18) 
(19) 



Calculating (cos S4>(t)) as J cos[A0 — A0o] 

<j{A(j), t\ A0 o ) cr st (A0 o ) d(A4>) d(A(f> ), we finally find the 
qubit correlation function 



v i s. COS^oT 

k z{t) = — exp 



C 



16F 



(20) 



The validity range of this result is the same as for Eq. 
iPfl (C < 1. 16F/C > 1); we have checked that in this 
range Eq. (12011 fits well the numerical Monte-Carlo re- 
sults. Fourier transform S z (u>) = 2 j ' _ oQ K z {f)e wl ' dr of 
Eq. (|2Ti|l in the case of efficient feedback (C/16F < 1, so 
the exponent is expanded up to the linear term) gives the 
oscillation spectrum (oj > 0) 



2tt 

i + f 2 + {uj/n ) 2 



_c 

~8Sh [TTf 2 - (uj/no) 2 } 2 + 4F 2 (u/n ) 2 



, (21) 



in which the first term (5-function) corresponds to syn- 
chronized non-decaying oscillations, while the second 
term describes fluctuations and for F <C 1 is peak-like 
near ui « Slo with the peak height of C/I6O0F 2 and half- 
width at half- height of F^Iq. [It is easy to check that 
/~5»(w)dw/27r = l/2.] 

Let us also calculate the correlation function of the 
detector current Ki(t) = ([I(t + t) - I ][I(t) - 7 ]). 
Following Ref. H3, we use Eq. J5J) to get Ki(t) = 
(AI/2) 2 K z (t) + K^t) + {AI/2)K z6 {t), where K € = 
(Si/2)5(t) is due to pure noise while the cross- 
correlation term K z ^(t) is due to quantum back-action, 
which shifts the phase <j> by — smcf>(AI / Si) £(i) dt as 
a result of noise £ acting during infinitesimal time dt 
[see Eq. I|13fl ]. Because of the feedback, the effect 
of this extra phase shift decreases (on average) with 



5 



time as 56(t) = - exp(-Ffi O T)[sin0(i)](AI/Sj) £(f) dt 
[see Eq. <|18[l ] an d the cross-correlation at t > 
can be calculated as K z ^(t) = (z(t + r)£(t)) = 
(cos[6(t) + tt T + 86{r) + <50(t)] Expanding co- 

sine up to the linear term in S6(t) [the linear expan- 
sion is the reason why it is sufficient to keep only aver- 
aged value S6(t) instead of the full distribution], we ob- 
tain K z s(t) = {£ 2 (t) dt)(AI/Si)exp(-FQ, T){sin[6(t) + 
n T + 54>(T)} sm[cj)(t)]), where (£ 2 (t)dt) = 5//2. Using 
symmetry of fluctuations leading to (sin^^r)) = (as 
above) and averaging over fast oscillations (sin [</>(£) + 
£Iqt] am[6(t)]) — (cosOo r )/2, we finally obtain 



AI 

K zd T ) = — (cosfioT)e 



-FfinT 



cos 86{t)). 



(22) 



Since expression for K z (r) has a similar structure [see Eq. 
Q17[l]. the corresponding terms of Kj(t) are combined to 
yield 



-FUot 



)K z (r), (23) 



where K z (t) is given by Eq. Ij20|l . 

To calculate the spectral density Si(ui) of the detec- 
tor current, we again expand the outer exponent of Eq. 
(|20|l up to the linear term [validity of Eq. (|23|l requires 
16F/C > 1]; then the Fourier transform gives 



Si(u) = Si 
Si C 



(A/) 2 



1 



C 



16F 

F 2 [l + F 2 + (u J /n ) 2 



2tt 



4 F [1 + F 2 — {u/n ) 2 ] 2 + 4F 2 (w/fi ) 2 



T4, (24) 



where the last (fourth) term T4 is the same as the pre- 
vious (third) term but with F replaced by 2F and with 
extra factor C/16F [actually, higher-order terms of the 
exponent expansion will lead to extra terms with F re- 
placed by 3F, AF, etc., and will slightly change the coef- 
ficients of the existing terms] . Notice that the 5-function 
in the second term of Eq. (|24|l is due to synchronized 
nondecaying oscillations, while the third term at F -C 1 
describes a peak with height (Si/8)(C/F) and half- width 
-Frio near u> ~ f?o- 

It is easy to check that the integral over all terms in 
Eq. fZty except pure noise Sj, gives the total variance 
of the detector current equal to (A/) 2 /4 [this also fol- 
lows directly from Eq. (|23[) ]. the same value as without 
the feedback^ Similarly to the non-feedback case, this 
variance would naively correspond to the qubit jump- 
ing between the two localized states, instead of oscillat- 
ing continuously [which would give twice smaller variance 
(A/) 2 /8]; in the Bayesian formalism this fact is under- 
stood as a consequence of non-classical cross-correlation 
between output noise and qubit evolution. 

Concluding this Section, let us emphasize again that in 
the ideal case the sufficiently strong feedback (16F/C ^> 
1) forces the qubit evolution to be arbitrarily close to 
the perfect coherent oscillations running for arbitrar- 
ily long time. In this case the feedback efficiency D 



approaches 100%, qubit correlation function becomes 
K z {t) — (cosf2oT)/2, in-phase quadrature component of 
the detector current becomes equal (A/)/4, and the cur- 
rent spectral density contains (besides the pure noise) 
the J-function peak at desired frequency Qq with vari- 
ance (A/) 2 /8, and also the narrow peak around (if 
C/16 <C F <C 1) corresponding to same variance (A/) 2 /8. 



IV. EFFECT OF IMPERFECT DETECTOR 
AND EXTRA DEPHASING 

Various nonidealities reduce the fidelity of the quan- 
tum feedback preventing D from approaching 100%. In 
this Section we consider the effects of imperfect quan- 
tum efficiency of the detector (77 < 1) and extra qubit 
dephasing with rate 7 en „ due to coupling to environment 
(see Fig. Both effects contribute to the total qubit 
dephasing rate 7 = ^ env + (i]^ 1 — 1)(AI) 2 /ASi in Eq. 
(0} and can be characterized by effective quantum effi- 
ciency of the qubit detection r\ e = [1 + A^Si / {AI) 2 ]^ 1 = 
\r)~ l + 4'f env Si/(AI) 2 }~ 1 or by effective relative dephas- 
ing d e = 1 /[{AI) 2 /ASi] =p-i-l + A lenv Si/{AI) 2 = 
r\~ x — 1; the physical meaning of d e is the ratio of qubit 
coupling to sources of pure (unrecoverable) dephasing 
and qubit coupling to the detector governed by the quan- 
tum (informational) back-action. 

Extra dephasing d e makes the qubit state non-pure; 
however, it is still perfectly monitored in a sense that 
p m (t) — p(t) (we assume that the magnitude of dephas- 
ing is known in the experiment and is used in the pro- 
cessor solving the quantum Bayesian equations; we also 
still assume e = 0.) Therefore, controller J7J with suffi- 
ciently large feedback factor F can reduce the phase dif- 
ference compared to the desired oscillations practically 
to zero. As a result, we would expect that the feedback 
efficiency D(F) should reach maximum (saturate) at in- 
finitely large F similarly to the ideal case shown in Fig. El 
however, this maximum will be less than unity. The sat- 
urating behavior of D(F) dependence is confirmed by nu- 
merical (Monte Carlo) calculations - see Fig.|3fa). Below 
we discuss the calculation of the saturated value D max 
at F = 00 [Fig.Efb)]. 

The evolution of a non-pure qubit state with Repi2 = 
(since e = 0) can be parameterized as pi\ — P22 — Pcosd, 
P12 = iP(sm 6) /2, where purity factor P is between and 
1. Using Bayesian equations © - ®> we derive evolution 
equations for P and 6 (in Stratonovich form): 



AI 



,AI 



P = — (1 - P z ) (— P cos 6 + £) cos 6 - 7P sin 2 6, (25) 



= 2 



Si 
Hfb 

h 



sin 6 AI , AI , . . sin 16 , , 

-p^(— + (26) 



Sufficiently strong feedback J7J) makes the phase 6 ar- 
bitrarily close to the desired phase 6 = VL$t (mod 2tt), 
so the feedback efficiency is practically equal to the av- 
eraged purity factor: D max = (P). To find (P) in the 
case of weak coupling C /f] e <C 1 , let us perform first the 
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FIG. 3: (a): Quantum feedback efficiency D as a function of 
feedback strength F for several values of quantum efficiency 
■q e of the detection, (b): Maximum feedback efficiency D max 
(at large F) as a function of r\ e . Dots show the Monte Carlo 
results for coupling C = 0.1, solid line corresponds to Eqs. 
1311 - 1321 . and dashed line shows the approximate formula 



averaging over oscillations and later the averaging over 
remaining slow fluctuations. It is easier to work with P 2 
than with P, so we start with evolution equation for P 2 
which is obtained from Eq. (J2SJ as dP 2 /dt = 2PP. It 
is easier to average P 2 over oscillation period using the 
Ito formffl because the noise £ causes correlated noise 
of (j>, and only in the Ito form the average effect of the 
noise is zero. Using the standard rule2Sil£ we translate 
the evolution equation for P 2 into Ito form: 



dP 2 (A/) 2 



1 - P 2 )(l - P 2 cos 2 <j>) - 2 7 P 2 sin 2 



dt 2Si 

+ (2A//S/)P(1 - P 2 )(cos0K 
then averaging over <p is trivial: 

^ = (A^ P 2 
dt 2St 1 A 2 



(27) 



(l-^ 2 )(l-V)-7^ 2 +- A/ " /1 " 2 



po-p*)s, 

(28) 

where is now a different white noise but with the 
same spectral density S$ = Si, so we do not change no- 
tation. 

A simple estimate of D max can be obtained from Eq. 
(|28|l by neglecting the noise term and finding stationary 
value for P, which gives^ 

D max « [1 + l/2r? e - V(l + V2^e) 2 - 2] 1 / 2 . (29) 



If we do not neglect the noise term in Eq. I|28|l . 
then P 2 fluctuates in time, and the stationary probabil- 
ity distribution a s t(P 2 ) can be found from the Fokker- 
Planck-Kolmogorov equation similar to Eq. 115fl (notice 
that varying diffusion coefficient comes inside the second 
derivative term) that leads to equation 

[ 7 P 2 - (1 - P 2 )(l - P 2 /2)(AI) 2 /2Si]a st 

+ [(A/) 2 /25 / ]^y[P 2 (l - P 2 fa st ] = 0, (30) 

which has analytic solution a st (P 2 ) = NG(P 2 ), where 



G(P 2 ) = (1 - P 2 



-5/2 



exp 



2(1 -P 2 ) 



(31) 



and N is the normalization factor. Stationary probability 
distribution for P can be found as a st {P) — 2Pa st (P 2 ), 
and calculating the average P gives us finally the feed- 
back efficiency 



Drnax — 



Jo P 2 G(P 2 ) dP 
J^PG(P 2 )dP' 



(32) 



Figure OJb) shows the dependence of the feedback ef- 
ficiency D max on the effective quantum efficiency of the 
detector r\ e = (1 + c? e ) _1 . Solid line shows the analyt- 
ical result (|32|l . dashed line shows approximate formula 
(|29|l . and the symbols show the numerical (Monte Carlo) 
results for D max (at sufficiently large F) for coupling 
C = 0.1. Notice that the lines for exact and approximate 
formulas are quite close to each other. 

Since at finite detection efficiency rj e the ensemble 
qubit dephasing is T — ?/ < r 1 (A/) 2 /4Sj, the weak coupling 
condition requires C/rj e ^1. As a result, the numerical 
results for C = 0.1 in Fig. |3| start to deviate (upwards) 
from the analytical result at r/ e < 0.03. For larger C the 
deviation starts even at larger rj e . Numerical calculations 
also show that at C/rj e > 3 the average purity factor (P) 
has a noticeable dependence on the feedback factor F 
((P) decreases with increase of F), while at C/rj e < 1 
this dependence is negligible. 

It is easy to see that in vicinity of the ideal case 
(r/e 1) Eq. I|29|l gives linear approximation D max w 
(1 + r] e )/2 pa 1 — d e /2 [exact solution (|32|) shows the same 
linear approximation]. This explains the corresponding 
numerical result of Ref. ^2- I n the opposite limiting case 
i] e <C 1, Eq. (|29|l is reduced to D max w \J2r\ e ; the exact 
solution (|32|l has a similar dependence but with slightly 
different prefactor: D max w 1.25y/rfe. Because of the 
square root dependence, feedback efficiency is still signif- 
icant even for large magnitude of qubit dephasing due to 
coupling with environment. For example, if coupling with 
dephasing environment is 10 times stronger than coupling 
with quantum-limited detector [d e = 10, r\ e = 1/11), 
then D max ~ 0.36, which is still a quite significant value 
for an experiment. 
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V. EFFECT OF e AND H DEVIATION 

In the ideal case we have assumed symmetric qubit 
(e = 0) and assumed that the exact value of tunneling 
parameter H is used in the processor. In this Section we 
analyze what happens if the qubit parameters e and H 
deviate from the "nominal" values e = and H = Ho as- 
sumed by an experimentalist and used in the processor. 
In this case the monitored value p m of the qubit den- 
sity matrix differs from the actual value p; and because 
of the mistake in qubit monitoring, the feedback perfor- 
mance should obviously worsen. [Both p(t) and p m (t) 
satisfy Eqs. with the same detector output I{t); 

however, "incorrect" parameters e = and Hq are used 
to calculate p m (t), while actual evolution p(t) is governed 
by actual parameter values s and H.] The desired evolu- 
tion is still pfi = (1 + cosf2oi)/2, pf 2 = i(sin£!o£)/2 with 
ilo = 2i?o / fi, which is used in calculation of feedback ef- 
ficiency D [notice that in the definition of efficiency i|10|) 
p d (t) is multiplied by the actual density matrix p(t), not 
the monitored value p m (t)\. The controller is still given 
by Eq. Q (we do not replace here H by Hq because this is 
more natural, for example, for control of the Cooper-pair- 
box qubit). Since the analytical analysis of the problem 
is quite complicated, in this Section we present only the 
numerical results of Monte Carlo simulations. 

Let us start with deviation of e (while H — Ho). Fig- 
ure^a) shows dependence D(F) for C = 0.3 and sev- 
eral values of e. One can see that for sufficiently large 
energy asymmetry e/H the feedback efficiency D is neg- 
ative at small F, while it is always positive at large F. 
For relatively small values of asymmetry (\e/H\ < 1) 
the dependence D(F) apparently saturates at large F, 
while at larger asymmetry [|e/-ff| ^ 1.5; not shown in 
Fig. E{a)] D(F) has maximum at finite F. [We can- 
not exclude the possibility that even for small e/H, 
D(F) also has maximum, but it occurs at too large F 
which cannot be analyzed by our code due to numer- 
ical problems.] Notice that the feedback efficiency is 
obviously insensitive to the sign of energy asymmetry: 
D(-e, H, C, F) = D(e, H, C, F). 

Solid lines in Fig. 0J D ) show dependence of D max- 
imized over F, on energy asymmetry e/H for several 
values of the coupling C — 0.1, 0.3, and 1. One can 
see that at small e/H the dependence D max (e/H) is 
parabolic (zero derivative at £ = 0), which means that a 
small energy asymmetry of the qubit decreases the feed- 
back efficiency very little. Zero derivative at e = is a 
natural consequence of the symmetry D ma x(—e/H) = 
D m ax(e/H) (because of this symmetry, we show only 
positive e/H). As seen in Fig. 0Jb), significant decrease 
of D max starts at smaller e/H for smaller coupling C. 
Rescaling of the horizontal axis by \[C makes the curves 
quite close to each other (see dashed lines in the Figure) ; 
however, we are not sure if the scaling D max {e/H\fC) is 
really exact at C — ► 0. 

The dotted lines in Fig. 01 b) show dependence 
D m ax{e / H) for a different situation, when the exact value 
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FIG. 4: (a): Dependence D(F) for several values of the qubit 
energy asymmetry e/H in the case when the processor and 
controller still assume e — 0. (b): Solid lines: maximized over 
F feedback efficiency D ma x as a function of the asymmetry 
e/H for three values of coupling C = 1, 0.3, and 0.1. Dashed 
fines: the same curves for C = 0.3 and 0.1 drawn as functions 
of e/HVC. Dotted lines: dependence D max (e / H) for the 
three values of C in the case when actual value of e is used 
in the processor, while the controller £J ls s tiU designed for 
e = 0. 

of e is used in the processor, but the controller is still 
given by Eq. @ designed for e — [desired evolution is 
still given by Eq. ijfjj with Sl = 2H/H]. One can see that 
exact monitoring of the qubit significantly improves the 
feedback efficiency compared with the case considered 
above; however, the feedback efficiency still decreases 
with energy asymmetry because the desired evolution © 
cannot be achieved at nonzero e/H and also because of 
non-optimal controller still designed for e = 0. (Some ap- 
parent dependence of the dotted lines on C even at C <C 1 
is possibly due to numerical problems of the code which 
does not work really well at C < 0.1.) 

To analyze the effect of the deviation of the parameter 
H from the value Hq used in the processor, we assume 
perfect energy symmetry, e = 0. Figure [SJa) shows the 
dependence D(F) for C = 0.3 and several values of the 
relative deviation (H — Hq)/H (we show only the curves 
for positive deviation; the curves for negative deviation 
are similar). We see that the effect of H deviation is 
qualitatively similar to the effect of energy asymmetry 
[compare Figs. EJa) and Eta)]. At large F the depen- 
dence D(F) saturates. FigureE^b) shows the value D max 
maximized over F as a function of the relative deviation 
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FIG. 5: Effect of the deviation of the qubit parameter H 
from the value Ho assumed in the processor, (a) : Dependence 
D(F) for several values of the relative deviation (H — Ho)/ H. 
(b): Solid lines: optimized over F feedback efficiency D max 
as function of the deviation (H — Ho)/ H for coupling C = 1, 
0.3, and 0.1. Dashed lines: the same curves for C = 0.3 and 
0.1 drawn as functions of (H — Ho)/ HC. 



(H — Hq) I H for several values of the coupling C. One can 
see that the dependence is almost symmetric for positive 
and negative deviation, and is parabolic at small devi- 
ation similar to the case of nonzero e discussed above. 
Also similar is the fact that weaker coupling C requires 
smaller deviation of H for the same value of feedback 
efficiency. However, the scaling with C is now different: 
the curves become close to each other if D max is plotted 
as a function of (H — Hq)/HC [see dashed lines in Fig. 
Etb)]. The different scaling is a natural consequence of 
the fact that small change of fi = V 'AH 2 + e 2 /h is lin- 
ear in H deviation but quadratic in e. The results pre- 
sented in Figs.0Jb) andOJb) can be crudely interpreted 
in the following way: D max decreases significantly when 
the Rabi frequency change due to parameter deviations 
(Afl = 2AH/H or AS1 ss e 2 /AHH) becomes comparable 
to the "measurement rate" (AI) 2 /ASj. [Notice that if the 
horizontal axis in Fig.[3Jb) was chosen as (H — Ho) /Ho, 
the curves would be somewhat more asymmetric, the 
asymmetry being more significant at larger C] 

Concluding the discussion of e and H deviations, let us 
mention that the main practical conclusion of the analy- 
sis is that the feedback operation is robust against small 
unknown deviations of the qubit parameters. 



VI. FEEDBACK CONTROL OF A QUBIT WITH 
ENERGY ASYMMETRY e 



In this Section we analyze the case of a qubit with fi- 
nite energy asymmetry e ( "asymmetric qubit" ) . In con- 
trast to the problem considered in the previous Section, 
in which nonzero e was treated as an unwanted deviation 
from the perfect zero value (therefore, finite e was just 
worsening the feedback designed for e — 0), now we try 
to design and analyze a different feedback (different con- 
troller) which goal is to maintain the free oscillations of 
a qubit with e ^ (so, now effect of nonzero e is what 
we also want to protect from decoherence) . Hence, the 
desired evolution p d (t) is no longer given by Eq. ©. 

Before choosing the desired evolution, let us mention 
that the qubit asymmetry leads to one more degree of 
freedom on the Bloch sphere. In case of e — 0, a pure 
qubit state was characterized only by the phase 4> [ see 
Eq. (|llfl ] because the real part of p\2 was vanishing in 
the course of measurement, so the evolution was within 
the plane of "zero longitude meridian" . For an asymmet- 
ric qubit, a naturally preferable plane of oscillations on 
the Bloch sphere no longer exists; in particular a weak 
measurement leads to a slow fluctuation of the "slanted" 
plane of free qubit oscillations (Fig.[SJ). The simulations 
show that without feedback the pure-state qubit evolu- 
tion is to some extent confined between the two slanted 
planes passing through the "north pole" (pn = 1) and 
"south pole" (p22 — 1), with the probability about 0.6 of 
being between the two planes for small C and \s/H\ < 1. 

Let us choose the desired qubit evolution as a free evo- 
lution starting from the north pole: 



P?i(*) 



2H 2 



e 2 + 2H 2 cos nt 



AH 2 + e 2 
= 1 + - cos 2 a (cos fit — 1), 



(33) 



pUt) = 



eH(cosQt - 1) 

AH 2 + e 2 
cos a 



iH sin Qt 



2 (AH 2 + e 2 y/ 2 

[sina(cosQ£ — 1) + isinfit]. (34) 



where fl = y/AH 2 + e 2 /% and a — atan(e/2iJ). An in- 
teresting question is whether the quantum feedback can 
keep the qubit evolution close to the desired path 
i|3l| or not. 

The old controller J7J) is obviously not good for this 
purpose, so we need to design a new one. [Notice that the 
qubit density matrix is monitored exactly, p m (t) = p(t), 
because all the parameters are assumed to be known ex- 
actly and because as discussed in Section II we assume 
infinite signal bandwidth.] As the first step, we charac- 
terize the deviation of the monitored qubit state p m (t) 
from the desired state p d (t) by two magnitudes (see Fig. 
|HJ): by the distance Ar m between two parallel slanted 
planes containing the monitored and desired states (the 
planes are slanted by angle —a) and by the angular dif- 
ference A(f> m between points p rn and p d projected onto 
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FIG. 6: Illustration of the qubit evolution on the Bloch 
sphere. For an asymmetric qubit (e 7^ 0) the free evolution is 
a rotation about a slanted (by angle a) axis. The difference 
between actual and desired qubit states (both are pure states) 
is characterized by the distance Ar m between the correspond- 
ing slanted planes and the angle Acj> m within the slanted plane 
(after projection). 



the slanted plane. The corresponding formulas are a little 
lengthy but straightforward. For the distance deviation 
Ar m = r m — r<i we calculate the distances r m and of 
the planes from the origin as the scalar products of the 
vector (cos a, 0, — sin a) orthogonal to the planes and the 
Bloch vectors (2Repi2, 2lmpi2, pn — P22) for the states 
p m and p d , correspondingly: 



2Rep™ 2 co s ol - (p™! - p™ 2 ) sin a 



(35) 



and similar for r^; it easy to see from Eqs. (I33(l - <|34[) that 
r c i = — sin a. For the phase difference A<fi m = <fi m — (fid 
(mod 2-7T, |A0 m | < 7r) we use equation 



tan < 



2I1M2 



2Rep^2 sin a + (p^ — p™ 2 ) cos a 



(36) 



[extra 7r-shift of <fi m is added when the denominator is 
negative as in Eq. ; a similar equation for p d defined 
by Eqs. - lj34"|) obviously gives (fid = fit (mod 27r). 
Notice that in the case e = (so that a = 0) we recover 
the previous definition © of A(fi m , while Ar m = 0. 

Limiting ourselves by the feedback control of the qubit 
parameter H only, we designed and analyzed the follow- 
ing controller: 



AHf b = -FHA(fi m - F r H sin cfi m Ar 



(37) 



The first term in this expression is the same as in the 
previous controller (JTJ and is supposed to reduce the in- 
plane phase difference A<fi m by changing the oscillation 
frequency, while the second term is supposed to reduce 
the inter-plane distance Ar m . The idea is that the change 
of Hfb = H + AHfb affects the angle of the slanted plane 
of oscillations, and when it is done periodically in phase 
with the oscillations (due to the factor sin0 m ) the inter- 
plane distance can be gradually reduced. 
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FIG. 7: Feedback efficiency for an energy-asymmetric qubit 
(e 7^ 0). (a): Dependence D(F) for several values of the 
ratio F r /F. Inset shows the same curves at small F. (b): 
Optimized over F feedback efficiency D max as a function of 
F r /F for several values of qubit asymmetry e/H and two 
values of coupling C . (c): Feedback efficiency D mal optimized 
over both F and F r , as a function of asymmetry e/H for two 
values of C. Dots show numerical results while the lines just 
connect the dots. 



Numerical calculations show that this idea works really 
well. Figure[3a) shows the dependence D(F) for several 
values of the ratio F r /F using as an example parameters 
e/H = 1, C = 0.3, and r\ e — 1. One can see that non- 
zero F r can significantly improve the feedback efficiency 
D. While at F r = the dependence D(F) saturates at 
large F, at non-zero F r the efficiency D has maximum at 
finite F. 

Figure shows the optimized over F efficiency 

D m ax as function of the ratio F r /F for couplings C = 0.3 
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and 0.1, and three values of energy asymmetry e/H. One 
can see that each curve has maximum at some value of 
F r /F. Notice also that at zero F r , the curves for different 
coupling C practically coincide, while at finite F r /F their 
behavior significantly depends on C, with larger D max at 
smaller coupling. 

In Fig. E{c) we show the feedback efficiency D max opti- 
mized over both F and F r as function of energy asymme- 
try e/H for two values of the coupling C. As we see, finite 
asymmetry e/H prevents efficiency D max from reaching 
100%. However, the difference 1 — D max decreases with 
decrease of coupling C, crudely proportional to C (except 
the region of small e/H, where the accuracy of our calcu- 
lations is possibly insufficient to distinguish the curves; 
unfortunately there is no simple way to estimate the cal- 
culation accuracy). Therefore, we guess that for any 
asymmetry e/H, the feedback efficiency D max reaches 
100% in the limit of small coupling C — » 0. (We cannot 
check this conjecture numerically because our code does 
not work well at C < 0.1.) 

Limiting the feedback control by the control of the pa- 
rameter H only, is not quite natural for the asymmetric 
qubit (though it is simpler from the experimental point 
of view). We have also performed a preliminary analy- 
sis of a simultaneous control of both H and e. We have 
considered the case when Eq. (I37f) is used for ii-feedback 
while a similar equation (with H replaced by e) is used 
for simultaneous control of e. Even though we have not 
performed detailed optimization, we have obtained larger 
values of D max than those presented in Figs. Efb) and 
dc). This shows that additional feedback control of the 
qubit parameter e really improves the feedback efficiency. 

Concluding this Section let us mention that its main 
result is the possibility of a very efficient feedback control 
of an asymmetric qubit. This can be done even using the 
control of the parameter H only, while simultaneous con- 
trol of e further improves the operation of the feedback. 



VII. CONCLUSION 

In this paper we have analyzed the quantum feed- 
back control of a single qubit, designed to maintain per- 
fect (or close to perfect) Rabi oscillations for arbitrarily 
long time. We have considered "Bayesian" feedback 12 
which requires a "processor" (see Fig.d solving quantum 
Bayesian equations to monitor the qubit state evolution 
via continuous output signal from the detector (QPC or 
SET) weakly coupled to the qubit. After comparing the 
randomly evolving (due to quantum back-action) moni- 
tored qubit state p m with the desired state p d , the qubit 
tunneling parameter H is being slightly changed in order 
to reduce the difference between the states. For simplic- 
ity we have assumed infinite bandwidth of the (noisy) 
detector signal and neglected the time delay in the feed- 
back loop. 

The analysis in Section III shows that in the ideal case 
the efficiency D of the quantum feedback can be made 



arbitrarily close to 100% by increasing the strength of 
the feedback control [characterized by parameter F in 
Eq. J7J]. It is important to mention that F scales with 
the coupling C between qubit and detector; therefore in 
the realistic case of weak coupling C -C 1, the parameter 
F and consequently the relative change of the qubit pa- 
rameter H remain small. The efficient operation of the 
feedback loop is achieved at F/C >■ 1 [see Eq. Q16fl]: in 
this case the qubit evolution becomes almost perfectly 
sinusoidal [Eqs. i(2"U)l and J2J], while the spectral density 
of the detector current [Eq. I|24fl ] contains the 5-function 
peak at Rabi frequency flo with the integral (A/) 2 /8 (as 
would be expected for the synchronized classical sinu- 
soidal oscillations in the qubit) and also a narrow peak 
around f^o with the same integral. The total integral un- 
der the peaks is thus (A/) 2 /4, which exceeds the limit 
for a classically interpretable process4i 

The feedback performance worsens in the case of a non- 
ideal detector and/or presence of dephasing environment. 
This case is considered in Section IV. We have obtained 
an analytical formula [Eqs. dH ] f° r ^ ne maximum 
feedback efficiency D max confirmed by Monte Carlo sim- 
ulations. It gives D max ps (1 + ?7e)/2 in almost perfect 
case when the effective detection efficiency r) e is close to 
unity, and D max w 1.25^/7^ when r] e < 1. 

In Section V we have analyzed numerically the de- 
crease of the feedback efficiency in the case when ac- 
tual qubit parameters e and H differ from the assumed 
(in the processor and controller) parameters e = and 
H = Ho (otherwise the case is ideal). We have found 
that for small deviations the efficiency D max decreases 
relatively slowly (with zero derivative at vanishing devi- 
ation), so that, for example, D max > 0.95 is possible for 
\e/H \ < 0.5VC and \H/H - 1| < 0.03C. This shows 
that the quantum feedback is robust against small devi- 
ations of the qubit parameters. 

In Section VI we have analyzed the feedback control 
of a qubit with finite energy asymmetry e, so that the 
desired evolution trajectory is along a slanted circle on 
the Bloch sphere. Despite the control problem becomes 
two-dimensional in this case even for a pure state, we 
have shown that efficient feedback is still possible using 
only one controlled parameter H and properly designed 
algorithm (controller) . 

The Bayesian quantum feedback of a solid-state qubit 
analyzed in this paper is not yet realizable experimentally 
at the present-day level of technology. Even much sim- 
pler quadrature-based quantum feedback^ is still a big 
experimental challenge. However, a rapid progress in ex- 
periments with solid-state qubit and also recent realiza- 
tion of quantum feedback in optics^ allow us to believe 
that the analysis performed in this paper will eventually 
be experimentally relevant. In this paper we have not 
considered two more effects quite important for the op- 
eration of the Bayesian quantum feedback: finite signal 
bandwidth and time delay in the loop. These effects will 
be a subject of a separate publication. 
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